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1 Introduction 

The calculation of fundamental solutions is difficult even for simple geometries. However, in the 
case of Laplace’s equation on the n-dimensional sphere, the sphere’s homogeneity allows us to 
reduce the partial differential equation to an ordinary one, rendering the computation tractable. 
The author is aware of two previous calculations of this type in the literature, [2] and [13], which 
disagree. In this paper, we shall proceed along similar lines to [2]; in the process, it will become 
apparent that the calculation in [2] is in fact the correct solution to a subtly different problem. 
We shall also check the consistency of our results with those of [13] in the case n = 2. 

The even- and odd-dimensional cases are treated separately: one may be expressed as a hyper¬ 
geometric polynomial in cot ^6, the other as a sum of an even polynomial in cotO and an odd 
polynomial multiplied by tt — 9. There is reason to expect this disparity, given the other ways 
in which spheres of even and odd dimension are different, such as the hairy ball theorem and 
the volume in terms of factorials. 

The last section gives some applications of our result, including the dipole potential on the 
sphere, the azimuthal Fourier expansion, fundamental solutions on real projective space by 
considering it as a quotient of the sphere, and a similar calculation explaining the problem 
that [2] is actually solving. 


2 The generalised Green’s function 

Let M be a compact manifold, S'” the (unit) n-dimensional hypersphere, and the hypersphere 
of radius R. 


*This paper is a contribution to the Special Issue on Orthogonal Polynomials, Special Functions and Applica¬ 
tions. The full collection is available at http://www.emis.de/jonrnals/SIGMA/OPSFA2015.html 
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In a recent paper [2] is what is claimed to be a calculation of a fundamental solution of 
Laplace’s equation on hyperspheres. The starting point of that paper is the equation 

—AG{x,x') = 6{x,x'), (2-1) 

where x,x' G and —A is the Laplace-Beltrami operator. However, this is not a well-posed 
equation on the sphere (or, indeed, any compact manifold)^; the best way to demonstrate this 
is to consider the Poisson equation that we normally use fundamental solutions to solve: 

—Art = p, (2.2) 

where p e n C^{S^). 

The divergence theorem. 


IM 


[diV X)dVg= / {X,N)gdVg 


IdM 


for any u G C°°{M) [10, Theorem 16.32]^ applied to gradn shows that 
[ {-A)u = 0, 


because dS"^ = 0; (2.2) then implies that 



If we interpret (2.1) distributionally, we should have 


{-A^{G{x;y),ip{y))y,x{x))x = {{G{x-,y), cp{y))y, -A^x{x))x 

= {{^xiy),^iy))y,x{x))x = {v’{x),xix))x = / ^x, 

Jm 

for every pair of test functions (p,x & ■= G°°{M), where the first equality is the definition of 

the Laplacian as a distributional derivative, the second the definition of G as a distribution, the 
third the definition of dx, and the fourth the definition of the pairing for smooth functions. 

If we take x(x) = 1 in the above equation, we arrive at a contradiction, since the first line 
must be zero and the last line need not be; hence, we need to modify the starting point, and 
consider a different initial definition of G. Suppose that we instead include a constant on the 
right-hand side of (2.1): 

—AG{x,x') = 6{x,x') — a. 


Carrying out the same procedure, we obtain 

{-Ax{G{x;y),ip{y))y,xix))x = {{dxiy) - a,ip{y))y,xix))x = / ^X-a{ ^]{ x] , 

JM \Jm / \Jm / 

which can only be consistent with the x(x) = 1 case if a = 1/|M|, where |M| is the volume 
of M, i.e.. 


|M| = [ dVg. 

Jm 

^On the other hand, the later paper [3] is unaffected, since the hyperbolic manifolds therein are not compact. 
^This is, of course, just a recasting of the fundamental theorem of exterior calculus, duj = 
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The behaviour is then that required of the usual Green’s function in the space of functions 
orthogonal to 1. 

See also [5, p. 354 ff.], wherein this is discussed as the notion of “Green’s function in the 
generalised sense”. Much of the theory carries straight over, provided that 5 is replaced by 
6 — |M|~^ in all cases^. 

Lastly in our general discussion, we note that a function G with these properties need not be 
unique; in particular, if we have two solutions, G, G', then 

-A(G - G') = (5 - \M \-^) - (<5 - \M\-^) = 0, 

so the difference is a harmonic function. The maximum principle implies that the only harmonic 
functions on a compact manifold are constant^, so G is determined up to an additive constant, 
which we may choose freely. Our convention shall be to choose the constant so that G is 
nonnegative on 5]^; we shall construct it to be so with a zero at the antipodal point to x. 


3 The ordinary differential equation for Gn 

3.1 The Laplace—Beltrami operator on a hyperspherical manifold 

To be explicit, let the Green’s function for the n-dimensional hypersphere be denoted Gn- One 
standard coordinate chart on the n-dimensional hypersphere is {6i, f), 0 < 6i < n, 0 < ip < 27r; 
this shall suffice for our purposes. (Our derivation does not depend on ip, and excluding the 
point antipodal to 0 = 0 will also not matter, since the function is constructed to be smooth 
there.) The metric is given by [9, p. 52, Example 7.4] 

g = diag sin^ sin^ 6 i sin^ 02, 


Taking 6i = 9, since a fundamental solution in a homogeneous set does not depend on 
direction, our fundamental solution reduces to G„ = G„(0), so we find that the Laplacian 
reduces to 

Af(9) = — 

^ ^ R^sin^-^ede^ ^ ’’ 

Therefore we need to solve the ordinary differential equation 

f + (n — 1) cot Of' = aR? , (3.1) 


for 9 G (0, tt), where a is the constant defined in the previous section. (Since the interval does 
not contain 0 = 0, the 8{9) term is not present.) The theory in the previous section dictates 
that, with Sn the volume of the unit n-sphere, 


1 _ r((n + l)/2) 

SnR^ 27r('^+l)/2i?« ’ 


(3.2) 


which is the reciprocal of the volume of the n-sphere; we discuss the relationship of this global 
constraint with the behaviour as 0 0 in Section 3.3 below. 

^See, for example, [1, p. 108 ff.]. 

'^See, for example, [16]. 
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3.2 Formal solution 

Recall that a first-order differential equation 
y +gy = u 

may be written in the form 
e~^{ye^y = u. 

It is then simple to show that for suitable limits on the integrals, 

y{x) = f dt 

Jxo 

is the solution, where xq is some constant. 

We conclude that (3.1), with a as in (3.2), has the formal solution 

GniO) = y / csc'^'V / yj dtp d(p. (3.3) 

‘JnU Je j<j) 

(Since the integrand is a smooth function of cp for 0 < (p ^ tt, the integral is smooth and bounded 
on [e,vr] for any sufficiently small e > 0 .) 


3.3 The behaviour as 0 4^ 0- consistency with fundamental solutions 
on flat space 


In this section we consider the limit for n 7 ^ 2; the two-dimensional case is analogous using 
log s. The integral (3.3) was determined by the global behaviour of fundamental solutions {viz. 
smoothness away from the pole at 0 = 0). On the other hand, the arc-length on the hypersphere 
of radius R is given by s = R9, so we expect to find that in the flat-space limit as i? —)• 00, 


Gn{s/R) ~ -r- 
[n 


1 

2)5„_is-2’ 


which shall give us the local behaviour as 0 i 0. We wish to verify that our solution has this 
local behaviour. To agree with the flat-space fundamental solution, our fundamental solution 
should also exhibit a singularity of order near 6 = 0, with coefficient 1/{{n — 2)Sn-iR^~‘^)■ 
We therefore need to have 


esc"" ^ (p I sin” ^ Ip dip d(p 


n—1 , 


SnR 


n—2 


(n-2)5„_i(R0)”-2 


as 0 i 0. We can drop the part of the (p integral from 7 r /2 to vr as irrelevant, it being just 
a constant, and then we need to examine 


f7r/2 


SnR 


n—2 


CSC 


„n-l 



Ip dip dip. 


It is apparent that as 0 j, 0, the integrand behaves as 


esc” p / sin” ipdp ^ p 
J<f) 


il—n 


sin” ^ p dp = p 


7rV2r(n/2) 

r((n + l)/2)' 


A result in asymptotic analysis shows that if F is sufficiently well-behaved and F{x) r\j fix) as 
X I xo, then for x G (xq, xq + e), e > 0, we have F ~ / as x | xq. It is easy to see 

that the integrand behaves as such, and hence we find 


R2-n 


f7r/2 


CSC 


n—1 



sin 


n-l 


p dp dp 


rsj 


02-nj^2-n 7 rV 2 r(n/ 2 ) 

(n-2)5„r((n + l)/ 2 ) 


(Rg)2-” 

(n - 2)5n-i ’ 


as expected, using the recurrence relation for the hyperspherical volumes. 
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4 Calculating explicit expressions for our fundamental solution 

In the following sections we shall consider two separate ways of determining our fundamental 
solution explicitly: the first proves that the solutions are finite sums of various trigonometric 
functions by deriving recurrence relations, the second uses hypergeometric functions and identi¬ 
ties to produce closed form expressions for these sums, and recurrence relations that determine 
the explicit form of the solution. Even and odd dimensions are treated separately since the form 
of the solution expression is quite different in each case. 

We shall use the notation 

:= [ sin™"^#, J^{0) := f , d(f). (4.1) 

J<j> Je sm™ ^ (j) 

It is apparent that Jm{d) shall be proportional to our fundamental solution: in particular, 

o2—n 

Gn{e) = ^Jn{ 0 ). 


4.1 Recurrence relation 

Proposition 4.1. The integral Im{4>) satisfies the recurrence relation 

Imifi) = ——H-^ sin™-”^ (/)cos (/) (4.2) 

m — 1 m — 1 

for m > 0, with basis cases 

= TT - (j), I 2 ((/>) = 1 -I- COS (/>. 

Proof. It is easy to verify the basis cases directly. We can carry out a standard calculation to 
find 


Imifi) = / sin*” — cos^ ff) dip = Im-2{4>) — / sin™' if cos^ if dip 


— Im—2{<f) 


1 


m — 2 


1 


sin™ ^ if cos if 


= Ira- 2 {(f)^ -^sin™ ‘^(fcoscf - ^-—Im{(f), 


+ 


1 


m — 2 


/ sin™“^ ■0(—sin^/:) dV’ 

J<t> 


I 


m — 2 


m — 2 


and we therefore have 


Im{(f) = —— jlm-2i<f) H-^ sin™ 2 . 

m — I m — 1 


>COS( 


for m > 0, as required. 

We may then verify, for example, 

h{(f) = ^-y^li(0) + ^ sin^“^ (fcosf = ^(tt - (f) + '^ sin 0 cos 0. 

We can use this to find a recurrence relation for Jm{d). 

Proposition 4.2. For m > 0, the integral Jm satisfies the recurrence relation 

1 cos 9 ^ I 

-2 n^'rn-2{fi) + 


'jyi _ Q 

Jm{e) = - -J^-2{e) + 

m — I 

with basis cases 


m — 1 sin™ ^ 


(m — I)(m — 2) 


J2i0) = logcsc^ ^6. 
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Proof. The basis cases are simple: for example, 


<^ 2 ( 6 *)= [ — dcj) = [ cot = logcsc^ 

Je sm (p Je ^ 

For the general case, [2] notes the useful trigonometric/partial fraction identity 


1 


1 

-.m— 2 


COS^ (p 
^ sin^d.’ 


sm " (p sm" 

which is a trivial consequence of Pythagoras’s identity. We therefore write 

Im{(p) cos^(p 


Jm{e) = 


le sm 


m-3 A ^ 


sm 


m— 1 


■ d(p. 


The next step is the most delicate part; we integrate the last integral by parts, being careful to 
note that near 9 = tt, Im{(p) ~ ~ (p)^- 


Id sm 


Im {(p) cos (p d(p 


1 Imi(p)cOS(p 
m — 2 sin”*“^ 1 


+ 


1 


1 ^( 0 ) cos (p - Im{(p) sin (p 


1 


- 


m-2 sin’"-^ 9 
— - "-^JmicP) - — 


m — 2 

T_/ 

m-2 Jg 


n 


sm 

(pcoscj) 


m—2 


d(p 


^m—2 


m-2 sin”^-^ 9 
We therefore have 
1 


m — 2 


sm 

sin cp cos (pd(p — 


dcp — 


1 


ImijpP) 


m — 2 Jg sm 

1 r ^rnicP) 


m — 2 Jg sin 


m—3 


„ d(p 

m-3^ 

d(p. 


Jm{6) = — ^({m-3)f d(p + . lm{0) - f sm(pcos(pd(p 

m — 2 \ Jg sm”^ <p sm™ ^9 Jg 


The first integral is substantially easier: it splits as 


Ig sm 


and hence 


^m{(p) , , m-2 

m-3 ^ d(P= - 


#+— 
m—3 ^ ^ 


sin™-2 


(pcoscp 


m — 1 Jg sm 

m-2 1 

:rJm—2{d) H - 
m — 1 m — 1 


m — 1 Jg sm 

/•TT 

sin (p cos (p dcp, 


m—3 


dcp 


Jm{e) = '^^Jm-2{0) + 


1 


,m-2 


m — 1 "" ' ' m — 2 sin™ ' {m — l){m — 2)Jg 

Computing the last integral, we obtain the recurrence relation 


sin p cos p dp. 


Yfi _ *3 

Jm{0) = - TJm-2{9) + 


1 


COS 0 ^ 

,m-2 a^rn[9) + 


1 


sm 


m — 1^'"' ' ' m — 2 sin™“^ 0 ' ' (m — 1) (m — 2) 

We can simplify the answer by using (4.2) to give 

T //ix m — ‘i^ 1 COS0 ^ 1 

Jm{0) — :.Jm-2{0) H 1 • m-2 J . , . ' 

m — 1 m — 1 sm 9 [m — l){m — 2) 

as required. ■ 

We shall now consider another way to obtain this recurrence, which will additionally produce 
closed-form results. 
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4.2 Hypergeometric functions 
4.2.1 n even: stereographic projection 

Note that stereographic projection M"" can be defined in terms of spherical coordinates 

{0,02, ■■■, 0n-i, on S'"’ and spherical coordinates (r, 6 * 2 , • ■ •, 0'n-i-, ^') on M" by 

r = cot ^0, 0[ = 0i, 2 ^ i ^ n — 1 , Lp' = Lp. 


The metric is transformed as 

d0‘^ + sin^ ^ , 


where is the metric on the (n — l)-dimensional sphere defined by { 02 , 03 , ■ ■ ■, 0n-i, v), 

and likewise {dYl^_i)‘^ in the primed variables. The Laplacian becomes 


A/ = 


1 


d 


(1 + r'^Y 

4r' 


sin’" ^ 0 80 

2\n Q 


^n—1 


+ 


1 


n—1 


sin^ 0 

df\ 1 


Ae„_i/ 


dr V (1 + ^ dr 


H— 


S' ,f, 

n — 1 


where As„_i is the Laplacian restricted to the (n — l)-dimensional sphere, and likewise ^ 
in the primed variables. 

Therefore in the new coordinates, the Green’s function is given by the integral expression 


r (1 + x^)"-^ / r j/" dy\ dx 
Sn Jo yJq (1 + y2)n y ) x ' 


(4.3) 


We observe that projection from the antipodal point is just r' = tan = 1/r. This fact is 
neatly captured in the r 1 —)• 1/r symmetry of the integrands in (4.3). Using this expression, we 
shall prove 


Theorem 4.3. The Green’s function for the n-dimensional sphere of radius R is 


Gn{0) = 


R 


2—n 


n — 2 9 ^ 

r^3F2 


Sn \n{n-l) 


1 , 1 , 2 -n/ 2 ^ 
2 , l + n /2 ’ 


—r 


—log(l + r^)^ , 
n — 1 / 


where r = cot ^ 0 . 

Proof. Consider first the inner integral in (4.3). This can easily be transformed to one of 
hypergeometric Euler type^ via the substitution u = (y/x)^. 


f 


dy 


= X 


u 


n/2-1 


/o (l + y2)"y JO 

We now use the Euler transformation® 


(1 + x^n) 


x^n) ^du = 


n 


2 F 1 


n, n /2 
_1 + n/ 2 ’ 


—X 


2 F 1 


a, b 


= ii-zy 


2 F 1 


c — a, c — b 


to write 


x(l + x2)” 


■n/ 2 , n 2 


ri, 1 - n /2 2I 

l + n/2' 

“11 2 2^1 

1 + X"' 

l + n /2 ’ _ 


®This integral was not actually discovered by Euler: the first unmistakeable linking of this integral to the 
hypergeometric series appears in the doctoral dissertation of P.O.C. Vorsselman de Heer [15, p. 10]. See also [6]. 
®[7, Sections 9-10 in particular]. 
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We can now use the contiguous relationship between F, F{a ), and F{b ) [ 8 , Section 7, 
equation (7)] 


{b — a){l — z)F = {c — a)F{a ) — {c — b)F{b ), 


(4.4) 


where F is shorthand for 2 Fi(a, b, c; z) and F(a~^) for 2 Fi(a + l,b,c;z) etc., to write this as 
a linear combination of a polynomial and x/(l + x^). Here, a = 1, 6 = 2 — n/2 and c = 1 + n/2, 
so 


(1 - n/ 2 )(l + x^) 2 F 1 


1 , 2 n/2 2 

I 1 + n /2 ’ 


= n /2 — (n — 1 ) 2 i^i 


1 , l-n /2 2 ' 

L l + n /2 ’ 


therefore 

X 1 


2 F 1 


n, 1 - n /2 

2 

n-2 

1 , 2 n /2 2 

\Fn/2 ’ 

— X 

- , .^X2Fi 

n[n — Ij 

, ; —X 

[ l + n/2 ' \ 


+ 


X 


2 (n - 1 ) 1 + a; 2 ' 


We can now integrate this directly to immediately obtain the general form of the Green’s func¬ 
tion, 


Gn(0)=4^ rf^LJ^X2F, 

Jo \n[n-l) 


n(n — 1 ) L 1 


1 , 2 n /2 _ 2 


— X 


+ 


1 


2 (n- 1 ) l+x‘^ 


dx 


= 4- 




2 —n 


Sn 

j^2-n 


n — 2 
2 n(n — 1 ) 


X 3 F 2 


1 , 1 , 2 -n /2 __ 2 

2, 1 +n/2 ’ ^ 


+ 


1 


n-2 2 TP 

-J 3F2 


Sn \n{n — 1) I 2 , l-|-n /2 


1 , 1 , 2 -n/ 2 , 2 


. 


4(n — 1) 
1 


log (1 -|- 


- x=0 


+ ^rrlog(l + r2) ) , 
n — 1 


as required. 


Corollary 4.4. Suppose n = 2m ^ 2 is a positive even integer. Then the Green’s function for 
the n-dimensional sphere of radius R is 


G2m{e) 


j^2—2m 

S2m 


( m — 1 n^ri, 1 , 2 — m 9 

~F) -^ 3^2 \-r 

\m[2m — 1) I 2, 1 + m 


+ 


1 

2m — 1 



(4.5) 


where r = cot ^6. 

In particular, notice that since 2—m is a negative integer, the hypergeometric term is explicitly 
a polynomial. 

To verify the local behaviour, we extract from (4.5) the most singular term as 1/r —)> 0. 
Firstly, we have r = cot ~ 2/0 as 0 —)• 0. As discussed in Section 3.3, the Green’s function 
should act like the fundamental solution for flat n-dimensional space near 1 /r = 0 , viz., 


1 


logs, 


n = 2 , 


(n 


1 

2)5„_i 


^2—n 


n 7 ^ 2 , 


where S'„_i = 2'k'^G jY{n/2J) is the volume of the (n — l)-dimensional unit hypersphere, and we 
recall that s is the arc-length, given by s = |x| in M”, and s = RO on hyperspheres S'^. 

For 52 (i.e., m = 1), we find that 


G 2 




which is consistent with Section 3.3. 
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For m ^ 2, the most singular term is the largest power of r in the hypergeometric polynomial, 
which has degree m — 2, and can be simplified as follows: 

R2-2m (l)^_2(lV_2(2-m)^-2(-l)"^r2("^-^)-^ 

S 2 m m{2m-l) ( 2 )m- 2 (l + (m - 2)1 

F(m + 1/2) 2r(m — l)F(m + 1) / r \ 2 m -2 
- 27r™+V2 r(2m + 1) \r) 

F(m — 1) mF(m)r(m + 1/2) 
y/n2mT{2m) 

^ F(m-l) 1 2^"^-^ ^ r(m) ^-2^ ^ 1 ^„2-n 

22”^ (ii6»)2™-2 {2m-2)2TT^ {n-2)Sn-l 

using the duplication formula for the F-function; this is also in agreement with Section 3.3. 


Remark 4.5. We can verify that some simple cases are in agreement with the corresponding 
results in [13]. Their notation differs somewhat from ours, and they have the opposite sign 
convention for the Laplacian (i.e., A instead of —A). Specifically, their even-dimensional Green’s 
function for the Laplace case (corresponding to L = 0 in their paper) is given in their paper’s 
equation (4.41) by 


^{2n-\-2) 


(n, n') 


1 


(2n -|- l)S'2n-|-2 
n! 


log 


1 — n ■ n' 


A (2n -2k- 1)!! • n') 

(2n + l)!! 52 n +2 ^ k{n-k)\ (1-n-n')^ + const, (.) 


where n is our m — 1, and 


—A, A -|- 2q: 1 — X 
a + l/2 ’ 2 

is the Gegenbauer function (of the first kind), which for A G Z is obviously a polynomial; we 
have also used that C'q“^(x) = 1 for a / 0 to simplify their expression into (4.6). It is easy to 
see that n ■ n' = cos 0 in our notation. Putting n = 0 in (4.6) gives the two-dimensional case. 




r(A -|- 2 q;) 
r(A + l)F(2a 


• 2.^1 



1 , 1 — COS 0 

^ log-^-h const. 


(4.7) 


whereas if we put m = 1 in our expression (4.5), we obtain 


G2{0) = ^Iog(l + cot2i0) = ^logcsci0 = -^log^^^, 


which agrees with (4.7) when we take the opposite sign (to agree with the difference in the 
Laplacians). 

Similarly, taking n = 1 in (4.6) gives the four-dimensional case, 



1 , 1 — cos 6 

^log^— 


1 1C 


( 1 / 2 ), 


COS 9) 




3^4 1 1 — cos 9 

1 — cos 9 cos 9 


+ const 




1 — cos ( 


1 — cos 9 
+ ^ cot^ \9 


+ const 


-|- const. 
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where the constants differ between the last two lines (essentially we have used 1—cos 9 = 2 sin^ 
again). Our corresponding result is 



which is once again in agreement when we take the opposite sign convention. 

4.2.2 n odd: a recurrence relation 

Unfortunately the simple hypergeometric polynomial formula found in the previous section does 
not extend to the odd case. Instead the Green’s function has a quite different character, as we 
shall see in Theorem 4.7. 

The discussion will be carried out for 7r/2 < 0 < vr initially, then the end result will be 
obtained by analytic continuation. 

Lemma 4.6 (reduction of J 2 m+i to hypergeometric form). Let 7r/2 < 0 < tt. With Jn{9) as 
defined in (4.1), we have 



where T = tan^ 0. 


Proof. We shall make the substitution 



2u^/2 tan <j){l + utan^ (/>) 


du 


in the definitions of l 2 m+i and J 2 m+i from (4.1). This also gives 


sin'i/j = tt^/^tan(/i(l + utan^ ^) 


Then, returning to the notation of Section 4.1 and taking a general odd n = 2m + 1, 


l2m+l (*/*) 


f 

J(j) 


'TT 


sin2”^t/)dV' = -^tan2™+i 



and the usual Euler integral produces 


l2m+l {fi) 



which can be converted to 


d-2m+l ifi) 



using an Euler transformation. 
We then have 
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Substituting 

tan^ (j) 


V = 


tan^ 9 
the integral becomes 

Jm{d) = 


= : T ^ tan^ 0, 


tan (f)d 4 ) = 


dv 




2{1 + Tv) 

1 , 1/2 


■-Tv 


dv. 


2{2m + l)jQ 1 + Tv [m + 3/2 
The same contiguous relation as before (see (4.4) above) allows us to write this as 


J2m+i{9) — 


T 

4m 


1 


1 


Fi 


1 + Tv 2m + 1 [m + 3/2 


1, 3/2 


■,-Tv 


dv, 


which integrates immediately to 

J2™+.(e) = /-(iog(i + r)- 


4m 


T 


2m + 1 


3 F 2 


1, 1, 3/2 
2, m + 3/2 


-T 


as required. 


This lemma provides a form of the Green’s function, but rather an unilluminating one: indeed, 
it is not even obvious that there is a finite form to the expansion. We shall now derive the 
following explicitly finite form. 

Theorem 4.7. Suppose thatn = 2m + l ^ 3 is a positive odd integer. Then the Green’s function 
for the sphere of dimension 2m + 1 is 


G2m+l{0) — 


where y = csc^ 6. 




*S' 2 m+l 2m 


Hence the solution can be expressed as an even polynomial in cot 9, added to an odd poly¬ 
nomial in cot 9 multiplied by tt — 9. 

Proof. We shall employ an argument using reduction formulae, based on the contiguous rela¬ 
tions for hypergeometric functions. 

We consider the hypergeometric expression in J 2 m+i derived in the previous lemma: write 


A{m) := -- 


3 F 2 


1, 1, 3/2 


-T 


2m + 1 [ 2 , m + 3 / 2 ’ 

By the 3 F 2 contiguous relation 

{b - a - 1)F = {b - l)F{b-) - aF(a+), 
where F = F{a, ... ; 6,... ■, z) [12, equation (15)], we know A satisfies 


A{m) — A{m — 1) = 


2T 


2 F 1 


1, 3/2 


;-T 


(2m + l)(2m-l) [m + 3/2’ 

Hence we can write down A{m) as a sum of hypergeometric functions: 

m—1 


A{m) = H(0) + ^ 


2T 


k=0 


{2k + 3){2k + l) 


2F1 


1, 3/2 _ 
A: + 5/2’ 


-T 


(4.8) 



























12 


R. Chapling 


and we also find that 

yl(0) = -^log(l + T), 

which cancels the other logarithm term, confirming our suspicion of the initial solution. This is 
still not obviously finitary. Now consider the term in the sum 

2T 

{2k + 2,){2k + l) 

We have the other contiguous relation for a hypergeometric function [8, equation (9)] 

(1 — a + (c — 6 — l)z)T = (c — a)F{a~) — (c — 1)(1 — z)F{c~), 
so after rearranging, B{k) satisfies 


1, 3/2 _ 
A: + 5/2’ 


9 2h — ^ 

B{k) = -■ ■ ■ ^ + y^^B{k - 1), 

^ ’ {2k + l){2k) ^ 2k ^ ^ 

where we have written y = 1 + 1/T = csc^(0). We also have 
arctan v/T \ 


R(0) = 21- 


Vt 


(4.9) 


Now split B{k) as B{k) = a{k)l3{k), where q;(0) = B{0), and a{k) solves the homogeneous 
version of (4.9): 

,,, 2A; + 1 ,, 

«(fc) = yo:{k - 1). 

We can show that 

«(fc) = /(^"^^/^)i?(o) 


(for example, by checking the recurrence relation and initial conditions). Now, /3 satisfies the 
recurrence relation 


2 2A; — 1 2 

= - {2k+l){2k) + = ~ (2t + l)(2fc) 


+ a{k)j3{k - 1), 


and hence we obtain the simple recurrence for (3'. 




2 

{2k + l)(2A:)a(A:) 


Then, since /3(0) = 1, we obtain 


m 


k 


i-E 

1=1 


2 

{21 + l){2l)a{l)' 


and since 


{2l + l){2l)(^ 
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we find that 


B{k) = a{k)-Y,: 


2a{k) 


^(2/ + l)(20a(0 

Now we can solve (4.8) to obtain 


= 2 


k — l/2'\ arctan-v/r^ 


k J 


r 1- 


Vt 


r, k (k-l/2\ 

k ) k-i 

2,1^ (i+i/2\ y ■ 
1=1 \ 1-1 ) 


Mm) = -- 


One can see from the expansion near 6 = tt that the correct branch of arctan (tan 9) to choose 
is 0 — vr, so we can substitute back to find that 

(«) = ^ g ( (‘' ~M) / (1 + (- «) cot.) - 1 g gg/-) . 


Hence, 


G2m+l{d) — 


j^l-2m 1 ™ 1 / A - l/2\ kn ^ 1 V- k^^) 

g(( k ji' (l + ('-«)“‘«)-3ggT^!' 


k-l' 


S2m+1 2,Tn 

as required. 

The most singular term in G2m+i(0) is 

j^l—2m 


1 fm-3/2\ 2 m- 2 . 1 1 /m-3/2 

^2S;( m-1 


^2171+1 2tti 

and then the coefficient of is 

r(m + l) r(m —l/2)7r 


1 


TT 


m 


2m-l ’ 


1 1 /m-3/2 

52m+i 2m V m - 1 


vr = 


27r™+i 2mr(m)r(l/2) 
r(m—1/2) r(m + l/2) 


47r™+i/2 (2m — l)27r”^+^/2 ((2m + 1) — 2 )S 2 m ’ 

as it should be to agree with the flat-space results. 

Remark 4.8. As in the even-dimensional case, we can compare this with the result in [13]. The 
three-dimensional case is given in that paper’s equation (4.23): 


— (3) / /\ ± I V ' i V . t, 

Gn {n,n) = arccos (—n • n ) -|- const. 

^ 47r2 Vl - n • n' 

Since arccos (—n • n') = vr — 0, substituting in gives 


n ■ n 


G^n\n, n') =- -(tt — 0) /,i + const =-^(vr — 0) cot 0 -|- const, 

47r2 I sin 01 47r^ 

since sin 0^0 for 0 ^ 0 ^ vr. 

Whereas, if we put m = 1 (and i? = 1) in our result we find that 


1 


cos 0 


1 


(4.10) 




1 


which agrees with (4.10) when the sign convention difference is understood. 

A similar calculation using [13, equation (4.34)] shows that the five-dimensional results also 
agree up to a constant, ours being 


G^{9) = cot 0 -F ^(5 -h 3csc^ 0)^ . 
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5 Applications 

5.1 The dipole potential on a sphere 

In Euclidean space, the conventional way to define a dipole potential from a fundamental solution 
is by considering two point charges with charge q and —q at a and P{d) a distance d apart. We 
then take the limit q —)• oo, d 0 , so that qd is constant, and the path has a well-defined tangent 
vector p when d becomes zero. To understand this on the sphere, it is advantageous to make 
this a little more formal; 


Definition 5.1. Let 7 ; (—e,e) —)■ M be a path on the manifold M, with 7 ( 0 ) = a and 7 ( 0 ) = 
p £ Ta{M). Let G{x,y) be the Green’s function for the Laplace-Beltrami operator on M. The 
dipole potential at a with dipole moment p is given by 


H{x, a,p) = lim 
t—>-0 


- G{x,a) 
t 


p ■ VaG{x,a). 


The limit definition indicates that this describes two point charges of charge being 

brought together; the second equality shows that H is well-defined by prescribing a and p. 
Moreover, since the total charge is always zero, the dipole potential is also consistent on compact 
manifolds, such as the hypersphere. Since two fundamental solutions are subtracted from one 
another, we also lose any dependence on the constant we took in the definition of our particular 
fundamental solution, so this is in a sense the dipole potential. 

We therefore need to find the gradient of G{x, a). Since this involves moving the base point 
of the coordinates, and hence different definitions of the angle 6, it is advantageous to start by 
recalling that on the sphere of radius R, 9 can actually be defined as a particular solution to 
the equation 


cos 6 = 


X ■ a 

~W' 


since it is the angle between the vectors x,a £ We then find that 

1 


Vn.9 = 


^X, 


R? sin 9 ’ 
and so the chain rule gives 


p ■ VaG{x, a) = GU0)ip • ^a)9 = • x), 

R^ sm 9 

and using the group of rotations on the hypersphere, we can always use the spherical symmetry 
to place a, p and x in the same three-dimensional space, with a along the axis 9 = 0 and p in 
the plane 9 = (p = 0, so that p ■ x = R\p\ sin 9 cos (p, and then 


p ■ VaG{x, a) = \p\ cos (p. 

We should now verify that this function has the correct behaviour as i? —)• 00 , in that it 
becomes a dipole at the origin in the Euclidean space with distance s = R9, 


^G'J9) \p\ cos(p 


\p\ cos (p 


S„R 


'n—1 


CSC 




(0) 


sm 


-^n—l 


9 ' d9' 


\p\cos(p 7 r^/^r(n/ 2 ) "■ |p| cos (p 

Sn{R9Y-^T{{n+ l)/2) ^ 
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similarly to the calculation in Section 3.3. It is easy to check that the dipole potential at the 
origin in M”’ is given by 


(p • V,) 


\x — a\ 


1 2 —n 


(n - 2)Sn-i 


I \—^ 1 —n 

\x\ p • X S 


IpI cos (f> 


a=0 


Sn-l 


Sn- 


n —1 


with which the limit agrees. 


5.2 Fourier expansion of our fundamental solution on 

The form of the n = 2 result means that we can quite easily derive an azimuthal Fourier 
expansion. We follow [4], with some simplifications. We have 

cosd{x,x') = cos 0 cos + sin 0 sin 0'cos (cj) — (p'), 

with x = (cos (/> sin 0, sin iji sin 0, cos 0) in standard spherical coordinates (and similarly x'), d{x,x') 
the spherical distance between them, and inserting this into our fundamental solution gives 

G 2 {x, x') = —— log 2 ^ ~ ^ sin0' cos {(p — cp')). 

We can find the expansion of this quite easily using the following: 

Lemma 5.2. Let A > B > 0 and t G M. Then 


2 (A — a? — cos kt (\/— B'^ — A 

- log (H + H cost) = log ^^ + ( - 


Proof. We begin with the well-known formula 


OO L j , 

9 , ^ COS kt 

— log (1 -I- — 2r cos t) = 2 ^ 

fc=i 


k 


k=l 


|r| < 1, 


B 


(5.1) 


which can by proven by applying Euler’s formula to the right-hand side and computing the sum 
directly. We then have 


2r \ / 2 \ cos kt 

- log ( 1 - ^ ^ cost) = log {l + A)+2y^ 

' k=\ 


1 _ 1 _ ^2 


k 


and it is easy to see that taking 


2r 


B 

A’ 


(5.2) 


— log ^ costj = logH — log (A + Boost). 

To find the right-hand side, we need to find the correct value of r in (5.2), but since |r| < 1, the 
only choice possible is 


1 _ 1 _ ^2 

the left-hand side becomes 
B 


r = 


\/H2 -B^-A 
B ’ 


and some brief algebra gives the result. 
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Theorem 5.3. Let G (0, 27r) and 6,9' G (0,7r). Then the azimuthal Fourier expansion for 
G 2 {x,x') in spherical coordinates is given by 


2'kG2{x, x') 


1 (1 —m)(l + M) ^ cos— (/>') /(I + m)(l — 

~2 4 ^ ^ k V(l + ^)(l-"i )) 


where M = max {cos 6 , cos 6 '} and m = min {cos 6 , cos 9'}. 

Proof. We can apply this directly, with 2A = 1 — cos 0 cos and 2B = — sin 0 sin 0'. Then 

— B‘^ — A |cos 9 — cos 9'\ — 1 + cos 9 cos 9' 

B sin 9 sin 9' ’ 

and squaring, square-rooting and applying trigonometric identities, this in turn becomes 




y (1 + M)(1 -m)’ 

where M is the larger of cos0, cos 9' and m the smaller. We hence obtain the result. 
We can write (5.3) most succinctly as 

2 TrG 2 ix, x') = log (esc Q9>) sec ( 56 »<)) -h ^ , 


k=l 


where 0> is the larger of 9, 9', and 0< the smaller. Taking the limit as i? —)> oo with s = R9 
constant gives 

°° cos k{4> — cj)') / R9^''^ 


2 TrG2{x, x') — log 2R —)• — log {R9>) + 


k=l 


k 


log s> -h ^ 


cos k{(l) — (j)') ( s< 


k=l 


k 


S> 


R9y 

k 


which agrees with the flat-space result quoted in [4], which we can also derive from (5.1)' 


5.3 A fundamental solution on 

The real projective space MP”" can be considered as the quotient of the sphere 5” by the group 
{1,-1}, where —1 is the antipodal map x i—)• —x, isomorphic to Z 2 '. MP"' = 5 ”/Z 2 . This is 
realised in our spherical coordinates by associating the points 

{9i,^p) ~ (vr - 9i,±-K + ip). 

Therefore our fundamental solution for S"" can be made into one for MP” by setting 
Gm>^{x,x') = ^{Gs"ix,x') + Gsr^i-x,x')), 
or in terms of angles, 

G]Rpn(0i) = -{Gs^{9i) + Gs"{t^ — 9i)). 

This clearly has all the correct properties of a fundamental solution, such as matching singular¬ 
ities and smoothness. 

^We had to subtract an extra constant due to our convention that G 2 is positive, which cannot happen in flat 
space as log is not bounded, above or below. 
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5.4 Interpretation of prior results 

In this section we shall discuss the results of [2], and in particular, how our construction enables 
a consistent explanation of some of the properties of their solution. We shall apply a similar 
procedure to the previous section, but consider the functions odd under the group action. 
Consider first the equation 

—AG(x,x') = 6 {x,x') — 6{—x,x'); (5-4) 


this is consistent on the sphere, because the integral of the right-hand side is zero. It is essentially 
the opposite limit of the two point charge case to the dipole considered above: the separation 
of the charges is now the largest possible. Linearity allows us to conclude that 

G{x,x') := G{x,x') — G{—x,x') (5-5) 

is a solution to (5.4), where G is the generalised Green’s function discussed in Section 2, since 

—AG{x,x) = [6{x,x') — {Sn)~^) — {d{—x,x') — = 6 {x,x') — 6{—x,x'). 

In fact, (5.4) is the equation that the Green’s function found in [2] satisfies: to see this we 
only have to observe its oddness under the transformation x i——x, or in our terms, 0 i—)• vr — 0. 
This is readily apparent if we consider the form [2, equation (29)] 


cos 9 2 Fi 


1/2, n/2_ 2, 


3/2 


; cos 9 


Hence it must have a singularity at 0 = vr of the same nature and opposite sign as it does 
at 0 = 0, where it was designed to have a singularity that produces a 5-function, so it must 
satisfy (5.4). 

An instructive and simple example is given by the nontrivial case in two dimensions: we have 

G{x, x') — G{—x, x') = — logcsc^ ~ ^ logcsc^ \ — 9) 

dvr ^ dvr ^ 

= — ^ logsin^ \9 + ^ logcos^ \9 = ^ logcot^ 
dvr ^ dvr ^ dvr ^ 

which is precisely the dimension-two Green’s function calculated in [2]. It also follows imme¬ 
diately that this satisfies (5.d). 

We now proceed to prove the relationship in the general case: we recall that our integral, 
from (d.l), is given by 


/•TT /•TT 

Jn{9) = / csc"'“^ 0 / sin”“^ 

Je Jd> 


^|J dip d4>, 


and Gobi’s is 


In{9) := 


dx 


I ’71 _1 ^ 

’g sm X 


up to a constant. 

Proposition 5.4. For 9 G (0,7r), we have 
Jn{9) - Mtt-9) = + 


(5.6) 
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Proof. We have 


Jn{0) = / CSC 

Je 


n—1 


sin” ^ ipdil) + I sin” ^ il) dip \ dcp, 


fvr/2 


n—1 , 


/7r/2 


so the difference is, using = /„ + 

rir—O 


JniO) - Jni’n-- 0) = / CSC 

Je 


n-1 


/7r/2 


1-77/2 \ 

sin”“^ Ip dip + / sin”“^ ip dip d(p, 


valid for 0 < 0 < tt when we recall the standard definition 

Ja 

in the bracket is the same as 


t/2 


Sin 


n—1 


Ip dip 


V^r(n/2) 
2 r((n + l)/ 2 ) 


/.“if b < a. The first integral 


using the Beta-function; the second is odd about cp = 7rj2. A simple way to see this is to change 
variables and notice that cos”“^ x is an even function, so it has an odd antiderivative. Since 
csc”“^ (p is even about (p = 7r/2 (again, a change of variables will show this easily), the second 
term’s integrand is odd overall, and hence the integral over a region symmetric about </> = 7r/2 
is zero. We are therefore left with 


^r(n/2) 
2r((n + l)/2) 


rn—O 

/ esc' 

Je 


<pd(p] 


this time, since the integrand is even, the answer is double the integral from 6 to vr/2, which 
gives the result. ■ 


With that resolved, we just have to check the constants. In our case, 

r ('fl'l - + f)/^) T 

Using (5.6), we have 


Gn{e) - - 0) = rf(l'+i(/2) (’^-(^) " 

7riGr(n/2) r((n + l)/2) r(n/2) 

r((n +l)/ 2 ) 27r(^+U/2i?n-2 " 2 W 2 R «-2 ”’ 

which is the same as Gobi’s solution. Hence the function in Gobi’s paper is related to ours in 
the way we specified in (5.5). 


A Characterisation of 2 F\ with half-integer coefficients 

We can use the contiguous relations for hypergeometric functions to derive a collection of rep¬ 
resentation results. 

Let n, ni,n 2 ,... be arbitrary positive integers, and /i,/ii,/i 2 ,... be arbitrary proper half¬ 
integers, i.e., h,hi G Z -|- 2 for all i. We also adopt an abbreviated notation for the ordinary 
hypergeometric function, 

a, b 

■,z . 

c 

Let also 1^{z) be the field of rational functions in z with integer coefficients, and V{fi{z)) be 
the Z(z)-vector space with basis {/^(z)}. Then we have: 


[a, b] c] = 2 F 1 
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Theorem A.l. There is the following classification 0 / 2 -Fi hypergeometric functions with half¬ 
integer and integer coefficients: 

1) if a G C is arbitrary, and c G C\{0, —1, —2,..., —n}, then [a, —n; c] G Q[^] is a polynomial; 

2) [ni, 71-2; ns] G V{1), i.e., a rational function of z; 

3) [/i,ni;n2] G F(l,^r^); 

4) [ni,n2;h] G F(1, (1 — arcsinz^/^); 

5) [hi,h2-,n]eV{lE{z),^K{z)); 

6 ) [/ii,n;/i2] G 14(1,2“^/^ arctanhz^/^)®; 

7) [hi, h 2 ] hfi\ G V (Vl — z, z~^l'^ arcsin z^!'^'), 

where K{z) and E{z) are the complete elliptic integrals of the first and second kinds respectively. 

This is a refinement and specialisation of Theorem 1.1 of [14] that is useful in the odd¬ 
dimensional case considered in Section 4.2.2; in our classification, the rational functions are over 
the integers since the parameters of the hypergeometric functions are rational. 

Remark A. 2. The Euler transformation shows that Cases 4 and 7 are the same. 

Proof. Case 1 is trivial. The other results follow from the contiguous relations and induction. 
These allow the expression of any one 2 F 1 in terms of two others with parameters integer steps 
away, with rational coefficients®. In particular, we may arrange circumstances so that 

F h’ 

2 -Cl 

C 

where {a,/3 ,7} C {1/2,1,3/2}. In particular, in the cases in the theorem, the forms of [a,j3,'fi\ 
are 

2 ) [ 1 , 1 , 1 ]; 

3) [1,1/2,!]; 

4) [1,1, 3/2]; 

5) [1/2,1/2,1]; 

6 ) [1,1/2 ,3/2]; 

7) [3/2,1/2, 3/2]. 

Evaluating the power series of the given hypergeometric functions with these parameters gives 
the result. ■ 
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®This notation is illogical, but common: see [11, p. 80]. 

®For more detail, see [14]. 


= P{z) 2 Fi 


a, j3 

5 

. 7 


+ Q{z) 2F1 


a — 1, fi 

5 

7 
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